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Abstract 

A unified approach valid for any wavenumber q, frequency u, and tempera- 
ture T is presented for uniform ideal quantum gases allowing for a comprehensive 
study of number density and particle-current density response functions. Exact 
analytical expressions are obtained for spectral functions in terms of polyloga- 
rithms. Also, particle-number and particle-current static susceptibilities are pre- 
sented which, for fugacity less than unity, additionally involve Kummer functions. 
The q— and T— dependent transverse-current static susceptibility is used to show 
explicitly that current correlations are of a long range in a Bose-condensed uni- 
form ideal gas but for bosons at T > T c and for Fermi and Boltzmann gases at all 
temperatures these correlations are of short range. Contact repulsive interactions 
for systems of neutral quantum particles are considered within the random phase 
approximation. The expressions for particle-number and transverse-current sus- 
ceptibilities are utilized to discuss the existence or nonexistence of superfluidity 
in the systems under consideration. 
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Degeneracy and Superfluidity. 
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1 Introduction 



Momentum correlations of infinite range are considered to be a fundamental property 
of a superfluid jH [2]. The complexity in evaluating the relevant finite-temperature 
expressions to investigate this property in any physical system has probably kept the 
suggestion almost dormant. For a one-component fluid, the knowledge of longitudinal 
and transverse particle-current static susceptibilities, x\\(q) an d X±{q), m the long- 
wavelength (q — > 0) limit is required to decide about the existence of long-range 
correlations. Since X\\{q) is exactly known for any wavenumber q from the famous 
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/-sum rule [3], the study of the property here boils down to the evaluation of x±(q) 
for q — > 0. We would like to investigate nontrivial systems for which the relevant 
susceptibilities can be determined in exact analytical forms. 

If one considers weakly interacting quantum systems and takes into account the 
interactions within the random phase approximation (RPA), the susceptibilities can 
be determined from knowledge of the corresponding expressions for a noninteracting 
system. Of course, for an ideal Fermi gas at T=0, the longitudinal response func- 
tions are well known, see e.g. [U Chap. 12]. For a Bose gas, Pines and Nozieres |3J 
Chap. 4.2, Vol.11] have discussed that x±(q)=0 at T=0, but it has been emphasized 
by Pitaevskii and Stringari [5j p. 97] in the context of an ideal Bose gas that evaluation 
of (T-dependent) x±(q) cannot be carried out so easily as that of X\\(q)- Density cor- 
relations at finite T in ideal Fermi and Bose gases have been discussed in Ref . [6] . It is 
worth mentioning, however, that although transverse-current correlations at finite T 
are available in semi-analytical and/or numerical forms for Fermi gases E] but their 
neat analytical forms are still lacking. 

The main purpose of this work is to present in a comprehensive manner a unified 
study of response functions of longitudinal and transverse particle-current density (as 
well as of number density) for all values of q, to and T for gases obeying Bose-Einstein 
(BE), Fermi-Dirac (FD), and Maxwell-Boltzmann (MB) statistics. Exact analytical 
forms are obtained for spectral functions as well as static susceptibilities. These dy- 
namical correlation functions are shown to be expressible in terms of polylogarithms; 
earlier investigations [9] and references therein] in terms of polylogarithms have been 
mainly for thermodynamic quantities. The expressions obtained by us would be use- 
ful in theoretical investigations of many physical properties enunciated in [8] and are 
utilized in the companion article [10], hereafter referred to as Paper II, to study the 
magnetic susceptibility in quantum gases of charged particles. Moreover, there have 
been interesting new developments in recent years [5] [Til E2j [13], [Hj in the study of 
low-density atomic quantum gases wherein the particles have short-range interactions 
and our results might turn out to be useful in situations where density and/or current 
correlations become accessible in experiments on such systems. 

The second purpose of this work is to utilize a q- and T-dependent expression for 
X±{q) to establish the relation lim^o X±(<?) < lim g _>o X\\(q) from which the existence 
or nonexistence of long-range correlations is deduced. It is found that while FD and 
MB gases possess correlations of finite range, and thereby constitute normal fluids, at 
all temperatures, the BE gas is a normal fluid for T > T c only. For T < T c , however, 
the BE gas has current correlations of infinite range even in the noninteracting case 
implying that the existence of long-range momentum correlations is only a necessary 
condition for a system to be superfluid. 

The outline of the paper is as follows. In Sec. [2] we introduce the definitions of 
various physical quantities required for our studies and briefly describe the procedure 
for their evaluation. In Sec. [3] exact analytical expressions for spectral functions 
are obtained and the results for the current response spectra are presented graphi- 
cally. The normalized current relaxation spectra are also plotted in this section. The 
particle-number and particle-current static susceptibilities are presented in Sec. H] 
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both numerically and analytically. These results are then applied in Sec. [5] to examine 
the issue of Bose condensation and superfluidity in noninteracting as well as interacting 
systems wherein the influence of interactions is considered within RPA. We summarize 
our findings in Sec. E] and some necessary details of our calculations are provided in 
the Appendix. 



2 Basic Considerations 



If the interactions between particles are negligible, the calculation of the number- 
density response function 



x(q,t) = —([N cl (t),Nl(0) 

with the number-density operator Nq= Ek<r a k o- a k+ q ,(j as well as of the particle-current 
response tensor 

X^(q,t) = ^-(KW,^'(0)t]) (2) 

with the particle-current density operator = Eka h/m (k a + q a /2) a k ^ak+q,^ (a, oc' - 
x, y, z) for a system of M identical particles contained in a box of volume V will sim- 
plify considerably. Here Nq(t) and J q (i) are the Heisenberg operators corresponding to 
A" q = Jd 3 r e~ iqr A r (r) and J q = Jd 3 r e _iq ' r J(r), respectively, a ko .(ak CT ) are the Fock-space 
operators which will create (annihilate) a particle in the state | kcr), and the angular 
bracket (...) denotes thermal average over the grand canonical ensemble, implying 



N. 



q =o 



--N. 



The many-particle Fock space averages in Eqs.flT]) and (j2j) will reduce to averages 
in the single-particle Hilbert space for a system of noninteracting particles. In par- 
ticular, for a uniform ideal quantum gas with Hamiltonian H= Ek<r £ k Qko- ^ an d 
^k = h 2 k 2 /(2m), it is straightforward to show that the above response functions of 
number density and particle-current densities will reduce to 



x{q,t) = 

which depends on g=|q| only, and 

2n 



2n 

— V C k sin 



Ak(q) 

h 



Xaaffat) 



A k (q) 



(3) 



(4) 



Here n=N /V is the overall number density, C^=J\fk/J\f denotes the thermal-average 
fraction of particles having momentum hk, vg(q) = (h/m) (k a + q a /2) and A k (q) = 
£k+ q - £k- Also, A^ k = E<x ( al a a ka ) and Ek Af k = A/". 



For the uniform (homogeneous and isotropic) systems studied here, 



q a q 



Xaa> (q, t) = X\\ (q, t) — — + x± (q, t) S, 



q a q a ' 



(5) 
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will have only two independent components, namely the longitudinal and the trans- 
verse. These components are given from Eq.Q in conjunction with Eq.Q as 

a, a' " " 

where the relation to x(<?> t) is reflecting particle-number conservation, and 



fy.rv' \ + 



5S*(")V^-( 



A k (q) 

h 



(7) 



with £=q • k/ (qk). 



Denoting by ip(q, t) any one of the response functions in Eqs.([3H7]), the corresponding 
dynamical susceptibility if)(q, z) for Qz > is defined as 

« (9 ,,)=i r ite ^( q , t) =r^^i ( 8) 

JO J-oo IT CO — Z 

with the spectral function (u real) given by 

] roo 



i>"(q,u) = 7: dt e^(q,t) = 3 lim^g, uj + ie) (9) 

Z J-00 <^0 

and for the related static susceptibility ip(q), one finds from Eq.fjHJ) 

tjj(q) = lim Mm Jj(q,uj + ie) = I — - , (10) 

e-»0cj — >0 J-00 7T W 

where the Cauchy principal value is to be evaluated. 

Applying Eqs. (|HJ) and (jUj) on the results given in Eqs.([3]) and (jlj), we obtain the 
dynamical susceptibilities 

x(q,z)=2nJ2C k fi(hz, A k (q)) , (11) 

k 

2 

Z 

X\\{q,z) = — + — x{q,z), (12) 
m g z 

nh 2 

X±(q, *) = "if £ C k #0- - £ 2 ) h fa A k (q)) , (13) 

and the corresponding spectral functions 

X \qM = 2n^C k / 2 (^, A k (q)), (14) 
k 

X|i(<?,^) = -2-x"(9jW), (15) 

II g 2 



nh 



Jd(q,uj) = —r £C k k\\ - i 2 ) h (fiu, A k (q)) 



m 



The functions f\ and /2 appearing in the above equations are 

A 



fi(s, A) 



A 2 



and 



7T 



f 2 (e, A) = - [S(e - A) - S(e + A)] . 



(16) 



(17) 



(18) 



Another quantity of physical interest, closely related to the current spectra, is the 
normalized Kubo relaxation function which determines the evolution of a small sinu- 
soidal initial current fluctuation j||,j_(r, t=0) = V|| ) j_cos(q • r) of the fluid according 
to 

j||,±(r,t) = (j>i±(q,t) V|| i _ L cos(q- r) . (19) 

Here the spectral function 0| ± (q, u) corresponding to <f>\\ t j_(g, t) is related to the current 
spectra by (see e.g. Ref.[2J Chap. 3.5]) 

-l 

(20) 



reflecting Kubo's identity. 



x",±(q, 



UJ 



-oc 
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The evaluation of the correlation functions requires knowledge of Ck given by 



n 9s 



A 



M eP £ >< - r]X ' 



(21) 



where g s —2s+l is the spin-degeneracy factor for particle's spin s, X=e /3fl is the fugacity, 
and i] = +1, —1,0 correspond respectively to a BE gas, an FD gas, and an MB gas. 
The normalization condition ^ k C\ — 1 yields the implicit equation 



C3/2 (v\)=nA 3 [l-C (T)} 



(22) 



that determines the chemical potential fi=fi v (n,T) of a uniform quantum gas. Here 

Cv{z) denotes the polylogarithm of order u, A= (27th 2 (3/ m^J ^ is the thermal de Broglie 
wavelength and 

C Q (T) = 5 V , x 6 (T c - T) [l - (T/T c ) 3/2 ] , (23) 

is the average fraction of particles occupying the ground state, wherein Sij denotes the 
Kronecker delta, Q(x) the unit step function, and T c = {2tiTi 2 /rnksj (n/g s ((3/2)) 2 ^ 3 the 
Bose condensation temperature; C(3/2) is the Riemann zeta function. 



3 Spectral functions 



The spectral functions may be calculated explicitly by making use of k-sum expressions 
derived in Appendix IA1 Choosing f(e)=2nf2(hu,e) in Eq. (1A.4D and 
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Figure 1: (Color online) Longitudinal-current response spectra of uni- 
form ideal quantum gases, Xl(<?*> T*)={2m / n-n) x'» (q, ^>) for q*=0A 
as a function of lu*: BE gas (blue), FD gas (red), and MB gas (dashed) 
at reduced temperatures T*=0.491(a), 0.164(b), 0.114(c), 0.104(d), 
0.055(e), 0.011(f) corresponding to T/T c =4.5, 1.5, 1.05, 0.95, 0.5, 0.1, 
respectively. 



f(e)=nh 2 /m 2 fzifou, s) which implies F{e)=nn% 2 / (2m 2 )[Q(e — huj) — Q(e + hu)] for the 
indefinite integral in Eq.l l A. 51) . the k-sums in Eqs.( fT4l) and (fIBT) are evaluated to give 
for the density response spectrum 

X(q,oo) = 2C (T)nnhoj5(h 2 u} 2 - efj +A Gi(x,b, rjX) , 

and for the transverse-current response spectrum 



A 

xKq, u) = — - G 2 (x, b, rjX) 
mp 



(24) 
(25) 



with x=(3hu/ (2y/pe^, b=JpV q /2 = qA/ (4^), 



G v (u,v,w) = — ku (e- {u - v)2 w) - C v (e- {u+v ^w) 
vw 1 v v 



(26) 



and 



.4 



g s X(3^ = [l-C (T)]ri\y/t 
4A 3 



n(3 , 



4Cs/2 (??A) 

wherein the rewriting of A follows from Eq. (1221) . and the recursion relation Cih-i( m 
fodz z _1 ( u (z) with Co (2) — z /0- — z )i see ; e -g- Ref- has been used. 



(27) 



-(x-bf 



one recovers from Eqs. 



Noting that G u (x,b,r]X) ^ (1/6) (l - e" 4 
and (125]) the known results for an ideal Boltzmann fluid (t>T=l/ V '(3m): 



nwir 



y2hqvT 



1 - e~ phw 
6 



exp 



{huj - e q f 
2v\h 2 q 2 



(28) 




Figure 2: (Color online) Transverse-current response spectra of uniform 
ideal quantum gases, Xt(9*' a; *5 T*)={2m / nir) Xj_(<7> w). Parameters as 
in FigfTJ 



^ (29) 

which will reduce to the corresponding expressions for a classical gas [15] , if the limit 
h — > is performed. 

For plotting purposes, we introduce the dimensionless temperature T*=k B T/e u , 
where £ u =/» 2 ^/2m and k u =2 (6ir 2 n/ gs) 1 ^ 3 denote system-dependent units of energy 
and wavenumber, respectively. We have T* = [6a/7t£(3/2)] 2 ^ 3 « 0.109 and, in the 
case of fermions with fcp as the Fermi wavenumber, k u = 2kp. We also introduce 
dimensionless reduced spectral functions, Xl,t( ( 1*^ T*) = {2m/n-n) x'n ±(q, uj), which 
depend only on the dimensionless parameters q*=q/k n , uj*=huj / 'e n and T*. 

The reduced spectra, which follow from Eqs.( fT2l . (|24l) and ( J25l) for uniform ideal 
quantum gases, are displayed in FigsJH and [2] as a function of uj* for one fixed q*=0A. 
Current response spectra are odd functions of frequency and hence they are displayed 
for to* > only. The graphs (a)-(f) in each figure refer to spectra calculated for 
decreasing temperatures ranging from high (T ^> T c ) to low (T T c ) values. It can 
be seen that due to T being sufficiently high in FigsiHa) and[2]^a), there is almost no 
difference between the bosonic and the fermionic spectral functions and, as expected, 
both nearly agree with the response spectrum obtained for the MB gas. In FigsJT^b) 
and (c), the Bose gas spectrum is seen to develop a precursor to the sharp free-boson 
excitation peak at u*=q* 2 =0.16, although the temperatures are well above T c . The 
spectral weight in the longitudinal boson-current spectra appears deformed due to 
thermal excitations which, owing to the o; 2 -factor in the expression corresponding to 
Eq. (1121) . get amplified more for larger uj. This mechanism, which is not effective in 
the case of the sharp 5-function resonance below T c reflecting boson excitations out 
of the condensate, gives rise to the peculiar separation of spectral weight visible in 
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to" 



to" 



to" 



Figure 3: (Color online) Reduced normalized longitudinal- 
current relaxation spectra of uniform ideal quantum gases, 
c/)1(q*,u*;T*)=2e u /(hTr)<p'/,(q,uj). Parameters as in FigHJ 



FigsUKcHe). 

As opposed to the longitudinal-current response, there is no contribution from the 
boson condensate to the transverse-current response (see Figj2]). This is due to the 
factor k 2 in Eq.( JT3l) which suppresses the (k=0)-term in the sum and leaves only 
contributions of thermally excited bosons with p=hh ^ 0. Hence, similar to the ther- 
mal contributions in case of the longitudinal-current response (FigHJ), the transverse- 
current response spectrum will vanish completely as T — > (FigJ2J), because all bosons 
condense into the zero-momentum state. The maximum value of ^±(g,u}), which is 
reached at uj=e q /Ti for low temperatures, can be seen to vanish as T 2 when temper- 
ature is decreased to zero. The interesting behavior found for the transverse-current 
response of the uniform ideal Bose gas has a noteworthy implication to be discussed in 
Sec. El 

We also calculated the normalized current relaxation spectra of longitudinal and 
transverse current fluctuations according to Eq. fl20|) . In Figsj3] and HI reduced spectra 
4 > l,t( ( 1*^ w*; T*) = 2e u /(hn) (j)" ± (q,cj) are displayed. As opposed to the response spec- 
tra, the relaxation spectra are even functions of u>*, which are normalized according to 
Cda/0£ iT (gW;T*)=l. 

The longitudinal spectrum 0£ for a BE gas (blue lines in Figj3]) develops a sharp peak, 
as the temperature is approaching T c from above, with an additional 5-function peak 
for T < T c . This peak can here be traced back to the same excitations of condensed 
bosons which are also responsible for the corresponding peaks in Fig{TJ According 
to Eq. (|19|) . these sharp peaks have the following physical implications. Any small 
longitudinal current fluctuation of large wavelength will perform undamped oscillations 
of very low frequency, i.e. the flow will sustain for long times. Near the surface 
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0.5 1 . 



(d) 



0.5 1. 



(b) 



0.5 1 . 



(e) 




0.5 1. 



(c) 



0.5 1 . 




0.5 1. 



Figure 4: (Color online) Reduced normalized transverse- 
current relaxation spectra of uniform ideal quantum gases, 
cp^(q* ,co*;T*)=2e u /(hTr) 4>'[(q,io). Parameters as in FigHJ 



of a fluid container, spontaneous current fluctuations across the surface will be of 
longitudinal nature resulting in the undamped creeping flow across the container surface 
as observed in superfluids. This relaxation behavior of longitudinal current fluctuations 
in a uniform ideal Bose gas below T c fits in with the implications of the transverse 
current response discussed in Sec. |5j 

Contrary to the longitudinal current relaxation, the transverse counterpart <fi^ be- 
haves surprisingly "normal". In fact, the transverse-current relaxation spectrum of a 
Bose gas is very similar, at all temperatures, to that of the Boltzmann gas (cf. FigH]), 
since the condensate does not contribute as mentioned earlier in this section. It is 
noteworthy to point out that at low temperatures the transverse relaxation spectra of 
BE and MB gases, as opposed to the FD gas, show a maximum at the same non-zero 
frequency, tu max ps e q /% (cf. Eq.fl2HD). 



4 Temperature— Dependent Static Susceptibilities 

We insert Eq.( TT5i) together with Eq.( l24l) into Eq.t jTO]) and then use dx (i(ae~ x2 ) = 
v^r (3/2(01) to demonstrate that 

x«(g) = - (so) 

m 

for all q and T, which reflects the famous /-sum rule (and which could be derived 
more easily by taking z — > in Eq. (fT2l) . of course). However, the evaluations of x(q) 
and x±(<l) are more tedious. The use of Eq.flTUJ) in conjunction with Eq. fflty and Eq. 
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(|25|) . respectively, gives 



X(q) = n(3 



C Q (T) 2A r°°dx„ . , 



2b 2 



n(3ir Jo x 



and 



X±(q) = -^5-f — G 2 (x, 6, 77A) 
nmp Jo x 



(31) 



(32) 



For A < 1, i.e. for negative chemical potential of the gas, the above integrals may be 
evaluated exactly as a power series involving error functions with imaginary arguments 
which we finally express as 



x(q) 



q^O 



n(3 
n/3 



inC (T) [1 - gb(r)] ^ (y\y f 3. 2 
^Co(T) 



(33) 



and 



x±(q) 



A 2 q 2 

f (iM) A 2 C i faA) 



n[l ~ Co(T)] ^ (yXf ( 3 
mC3(r/A) £f /I 11 V 2' 



(34) 



£=1 



^ -[1-C (T)]<!1 



A 2 Ci foA) 
24tt (77A) 



g 2 + 0(q 4 



(35) 
(36) 



with iF% (atx, «2; 2) as the Kummer function, one of the confluent hypergeometric func- 
tions [15]. For an MB gas (r/=0), only the first term (£ = 1) of the series survives 
leading to closed analytical forms for both x(q) and x±(q). 

Whereas for a Bose gas the above expressions are valid at all temperatures, for a 
Fermi gas numerical evaluation of Eqs. (I3T1) and (13"2"1) is required in the regime A > 1, in 
general. We note, however, that the asymptotes given in Eqs. f[3"4"l) and (l3"6l) reproduce 
the small-g behavior of the susceptibilities correctly not only for BE and MB gases 
but also for a FD gas at all temperatures. The asymptotic Eqs. (134"]) and (136]) were 
deduced with the help of the Kummer-f unction series expansion 1F1 (l;3/2; — £b 2 ) = 
J2^Lo (—£A. 2 q 2 /8ii) p / (2p + 1)!!. But a word of caution needs to be in place regarding 
any temptation to use the latter in combination with a subsequent change of order of 
summations in Eqs. (13U]) and (133]) ; such a procedure would result in a divergent power 
series in q. 

It is seen that contrary to X\\(q) m Eq. (130]) . x±(q) depends on both q and T. Figure 
[5] displays the ratio XT—X±(q) / X\\(q) as a function of Q=T/T C for a set of wavenumbers 
in the range < q* < 1 for all the three gases, and with T c as the value corresponding 
to a uniform ideal Bose gas. Figures 0(a) and (c) show that the BE and MB gases 
have rather similar behavior for T — > except that the former shows a thermodynamic 
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phase transition at a finite temperature. The values X±(q) 7^ for a Fermi gas at 
T = in Fig. [5] (b) can be traced back to the fact that all states fill completely up to 
the Fermi energy. As T increases, x±{q 7^ 0) increases and asymptotically approaches 
X\\{o) — n/m irrespective of the statistics signifying that each of the gases considered 
here responds identically to transverse and longitudinal probes at high temperatures. 



Since lim^o G„(u, v , w)=Au C„_i [we " J /w, one finds from Eqs.(l30l and ( 1321) that 
uniform ideal FD and MB gases at all temperatures, and a BE gas at T > T c , will obey 

Ti 

lim-X±(g) = hmX||(g) = — ■ (37) 

q^O q^O m 

On the other hand, for bosons at T < T c we have 

11 

liiWg) = -[1 - C (T)] < limxn(g) , (38) 
q^o m q^o 

with the finite condensed fraction, Cq(T) > given by Eq.([23l. The inequality in 
Eq. ( 1381) achieved through explicit evaluation of transverse and longitudinal particle- 
current susceptibilities generalizes the T=0 result mentioned in Ref. [3, Chap. 4, Vol.11]. 



5 Condensation and Superfluidity 



This section constitutes a spin-off of our investigations in earlier sections. We wish to 
apply Martin's criterion to uniform quantum gases and closely follow the elaborations 
given in [2]. The /-sum rule given by Eq. fl30l) is a direct consequence of the gauge 
invariance and hence it holds for any fluid. The total mass density p = mn of a 
one-component fluid, interacting or noninteracting, is therefore defined as 



, f 00 du xf(?,o>) 
p = m x\\ [<1 — > 0) = m hm-h . (39) 



One defines another mass density p n by the relation 



Pn = m 2 x± (q -> = m 2 lim I — A±W ' ; 40 

and one finds p n = p from Eq. (1371) . implying momentum correlations of finite range 
which is characteristic of a normal fluid. However, the result from Eq. (!38|) . p n < p, 
implies that the gas has momentum correlations of infinite range which is a character- 
istic of a superfluid [TJ Sec. V.E],[2J Chap. 10]. Thus the difference p s =p — p n , known 
as superfluid density, turns out to be equal to the condensate density pC (T), a result 
also reported in Ref . [17] . 

We provide an independent consistency check for the expression p s = pCo(T) using 
the Josephson relation [18] rewritten in the form 



- dco G"(k,u) _ {m/hf 
I'Kfi oj p s k z 
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Figure 5: (Color online) Reduced transverse-current susceptibil- 
ity XT = X±{Q) m / n as a function of the temperature ratio 0=T/T c 
for uniform ideal quantum gases obeying BE (a), FD (b), and 
MB (c) statistics. In each plot, the curves refer to wavenumbers 
q*=0.0, 0.125, 0.25, 0.5, 1.0, respectively, from the top dash-dotted line 
to the bottom line. 
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which relates the superfluid (mass) density p s to the condensate (number) density 
n =nCo(T). Here G"(k,u) denotes the spectral function 



G (r - r', * - t') :=-£( [a CT (r, t) , 4(r', t')} ) 



(42) 



corresponding to the correlation function of boson field operators. For a uniform ideal 
Bose gas, one finds G"(k, uj)=ii5(uj — and hence obtains from Eq. (j4"Tl) . p s =mn 

in agreement with our finding. 

We note that the existence of long-range correlations together with our discussions in 
Sec. |3]of long-wavelength longitudinal-current relaxation spectra lead to the conclusion 
that a uniform ideal Bose gas at T < T c possesses some characteristic properties of a 
superfluid. However, there is a well-known result, which can also be retrieved from 
Eq. (1531) . that the static compressibility [15] given by K=lim g _ >0 /t(g) = lim g _» [x(q)/n 2 ] 
diverges for T < T c for an ideal Bose gas implying a sound velocity c=0 in the condensed 
phase. Hence such a system does not satisfy Landau's requirement for superfluidity, 
namely that of a nonzero critical velocity (see e.g. [5], EH]). This in turn suggests that 
the existence of long-range momentum correlations will only represent a necessary 
condition for a system to be superfluid. On the other hand, an interacting Bose gas, 
howsoever weak the interaction might be, has finite compressibility and sound velocity, 
and would be capable of sustaining a superflow. 

We, therefore, consider now particles interacting repulsively via a short-range pair 
potential v (|r|) whose influence can be incorporated approximately on our findings of 
previous sections by applying the RPA. The expressions for the dynamical susceptibil- 
ities given by Eqs. (fTTj) to ffTBT) then take the forms 



While Eq. (1431) is a well known RPA result [3], vol. I] and Eq. fT44|) follows from it as an 
immediate consequence of number conservation (cf. Eq. f[T2l) ). Eq. fT4"5"j) has been de- 
duced from Eq. (12) of Ref. [2D]- We note that Eq. fT4"5l) can serve as an approximate 
expression for the transverse particle-current susceptibility of an interacting system 
of neutral quantum particles only. For charged particles, transverse electromagnetic 
shielding effects will modify Eq. (1451) due to intimate coupling between particle current 
and charge current, as discussed in Paper II [10]. It may be remarked that the con- 
sidered approximation is a consistent weak-coupling theory having its validity in the 
limiting situation v(q) — > for short-range interactions. 

Assuming the pair potential to be a contact interaction, one has v(q)=4:iih 2 a/m, with 
a > denoting the scattering length. The compressibility of the interacting Bose gas 
turns out to be finite in the condensed phase, which is readily seen by taking the static 
limit in Eq. (l4~3~]) and using x{q) from Eq. (1311) . This removes the pathological divergence 




(43) 



(44) 



(45) 
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of the compressibility of the non-interacting gas and results in a finite sound velocity 
c, a necessary condition for the superfluidity of the Bose-condensed phase. It follows 
from Eq. (1441) that although RPA modifies the longitudinal dynamical susceptibility, 
its static value remains unaltered, as expected. Furthermore, Eq. (I45p implies that, 
within the RPA, interactions do not alter the dynamic, and hence the static, transverse- 
current susceptibility. This, therefore, implies that the relation between superfluid 
density and condensate density remains the same as in the ideal Bose gas. 



6 Summary 

We have presented a comprehensive study of the dynamics of uniform ideal quantum 
gases in terms of the number-density and the particle-current density response func- 
tions. Unified expressions for the associated response spectra, valid at any wavenumber 
q, frequency u, and temperature T, as well as the corresponding static susceptibilities 
have been obtained for BE, FD and MB gases. All results, many of them given in 
exact analytical forms and believed to be new, have been discussed in view of their ap- 
plications in various contexts. The plots of the response spectra provide an immediate 
comparison of the dynamical behaviour, indicating intrinsic and apparent differences 
due to the statistics. The relaxation spectra of longitudinal particle-current fluctua- 
tions have been found to develop a peak as T is approaching T c from above and lowered 
below T c , contrary to their transverse counterparts. The transverse-current relaxation 
spectra of the BE gas are surprisingly similar to those of an MB gas for all T. 

It is found in the long-wavelength limit that the static transverse boson-current sus- 
ceptibility approaches a T-dependent value equal to p n /m 2 < p/m 2 below T c . From 
this result we conclude that in a condensed ideal Bose gas the momentum correlations 
are of long range and the superfluid density equals the condensate density. Repulsive 
contact interactions have been taken into account within the RPA. This approximation 
removes the divergence of the compressibility and leads to a transverse particle-current 
susceptibility equal to that of an ideal gas. Our studies demonstrate that the existence 
of long-range momentum correlations constitutes only a necessary condition for a sys- 
tem to be a superfluid. 
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Appendix 



A Evaluation of k-Sums 



We have to evaluate the k-sums appearing in Sec. [2] of the forms S = yj k Ck/(A k (q)) 
and Sj_= yj k Ck k 2 {l — £ 2 )/(A k (q)) . Separating out the term k=0, converting the 
remaining sum over into the k-integration corresponding to a macroscopic system, 
and transforming to integration variable x=y/j3Fk, we obtain 

S = C„m/( £ ,) + J^jy x ^ r J\f(e q + 2^) , (A.1) 

and 

For later convenience, ^-integrations have been extended over the complete real axis 
since both integrands are even functions of x. Using the identity 

e x2 — r/X 2r]X x dx ( ^ ' ' ^ ^ 

the x-integrations are done by parts before evaluation of the ^-integrals is attempted. 
This allows Eqs. (lA.lj) and (1A.2I) to be rewritten finally as one-dimensional integrals, 

S = C (T) f(e q ) - J°°Jx In (l - V Xe' x2 ) f [e q + 2x^JfP) (A.4) 

and 

S± = xbx(l- r,Xe^) F (e q + 2xJ7jp) , (A.5) 

where F(x) = Jdx f(x) is the indefinite integral of the function /. 
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